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We present majority and minority radiofrequency (RF) spectra of strongly interacting imbalanced 
Fermi gases of 6 Li. We observed a smooth evolution in the nature of pairing correlations from pairing 
in the superfluid region to polaron binding in the highly polarized normal region. The imbalance 
induces quasiparticles in the superfluid region even at very low temperature. This leads to a local 
bimodal spectral response, which allows us to determine the superfluid gap A and the Hartree energy 
U. 



PACS numbers: 05.30.Fk 03.75.Ss 32.30.Bv 67.90,+z 

Pairing and superfluidity in fermionic systems are intri- 
cately related phenomena. In BCS theory [1 , describing 
conventional superconductors, the emergence of super- 
fluidity is accompanied by the opening of a gap in the 
excitation spectrum of the superfluid. This gap can be 
interpreted as the minimum energy required to break a 
Cooper pair or, equivalently, to create an elementary ex- 
citation, a so-called quasiparticle, inside the superfluid. 

However, strongly correlated systems show a more 
complicated behavior. There are gapless fermionic su- 
perfluid systems, as is the case for high temperature su- 
perconductors [2] or for superconductors with magnetic 
impurities [3 . On the other hand, there are numerous ex- 
amples of systems with an excitation gap in the normal 
state, e.g. a high temperature superconductor above its 
superfluid transition temperature exhibiting a pseudogap 
[2] or a semiconductor. 

Here, we use radiofrequency spectroscopy to investi- 
gate the nature of pairing and the relation between pair- 
ing and superfluidity in a strongly interacting system of 
ultracold atomic Fermions. 

We can spectroscopically distinguish the superfluid 
and the polarized normal fluid by introducing excess 
fermions into the system. In a superfluid phase de- 
scribed by BCS theory, the excess particles can be ac- 
commodated only as thermally excited quasiparticles. A 
double-peaked spectrum reflects the co-existence of pairs 
and unpaired particles. In the normal phase, at large 
spin polarization, the limit of a single minority particle 
immersed into a Fermi sea is approached, which can be 
identified as a polaron [4] [5j [6j [7] . Here the system can 
be described in the framework of Fermi liquid theory and 
no stable pairs exist. We find that these different kinds 
of pairing correlations are smoothly connected across the 
critical density imbalance [8], also called the Clogston- 
Chandrasekhar limit of superfluidity [9j [10] . 

Excess fermions in a low-temperature superfluid con- 
stitute quasiparticles populating the minimum of the dis- 
persion curve. The RF spectrum of a superfluid with 
such quasiparticles shows two peaks, which, in the BCS 



limit, would be split by A, the superfluid gap. There- 
fore, RF spectroscopy of quasiparticles is a direct way to 
observe the superfluid gap A in close analogy with tun- 
neling experiments in superconductors [11]. From the 
observed spectrum we can also determine a Hartree term 
P2], whose inclusion turned out to be crucial. 

For this study, we have combined several recently de- 
veloped experimental techniques: The realization of su- 
perfluidity with population imbalance [13] leading to 
phase separation [SJ [13j [14], tomographic RF spec- 
troscopy [15] , in-situ phase contrast imaging with 3D re- 
construction of the density distributions [8]. In order to 
minimize final state effects [16] we have prepared an im- 
balanced mixture of states |1) and |3) of 6 Li (correspond- 
ing to \F = 1/2, m F = 1/2) and \F = 3/2, m F = -3/2) 
at low field) in an optical dipole trap at a magnetic field 
of B = 690 G, at which there is a Feshbach scatter- 
ing resonance between the states |1) and |3) [16] ITT] . 
Evaporative cooling at B = 730 G is performed by low- 
ering the power of the trapping light. After equilibra- 
tion an RF pulse was applied for 200 /is selectively driv- 
ing a hyperflne transition from state |1) or |3) to state 
(|2)|F = 1/2, m F = -1/2) at low field). The RF power 
was kept constant for all experimental results presented. 
Immediately after the RF pulse an absorption image was 
taken of the atoms transferred into state |2). 

The spectra were correlated to the local Fermi energy 
e_p| = (67rn|) 2 / 3 of the majority density and to 
the local polarization g\ oc = which is a measure 

of the local excess fermion population. As in a previous 
publication [8] the local densities were measured using 
phase contrast imaging and 3D reconstruction using the 
inverse Abel transformation. 

The RF spectra shown in fig. [T] reveal a gradual 
change in the nature of the pairing correlations. The 
balanced superfluid is characterized by identical spectral 
responses of majority and minority particles and has been 
the subject of previous studies, see p~8j [19] [20] and ref- 
erences therein. In the polarized superfluid region [5J I2T] 
(and references therein) the minority spectrum perfectly 
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FIG. 1: (Color online) Tomographically reconstructed RF spectra for various regions of the atomic sample at unitarity. a) 
Balanced superfluid, b) polarized superfluid, c) moderately polarized transition region and d) highly polarized normal region. 
The panel on the left shows a phase contrast image of the atomic cloud before RF excitation. The positions of the spectra a) to 
d) are marked in the phase contrast image and by the arrows in figure [2] which displays all our results in the unitary limit. Red: 
Majority spectrum, blue: Minority spectrum. Local polarizations <ji oc and local temperature T/Tp, respectively: a) -0.04(2), 
0.05(1), b) 0.03(1), 0.06(l)c) 0.19(1), 0.06(2), d) 0.64(4), 0.10(2). The negative value in a) implies that the local polarization 
as inferred from phase-contrast imaging underestimates ai oc by up to 0.05. 




matches the pairing peak of the majority spectrum, lo- 
cally coexisting with the quasiparticle spectral contribu- 
tion, resulting in a local double peak structure of the 
majority spectrum (see fig. [T]d) . The spectrum suggests 
that the majority population can be divided into two 
distinct parts: One part consisting of pairs forming the 
superfluid, the other part consisting of quasiparticle ex- 
citations in the form of excess fermions. Therefore, a 
natural interpretation of the RF spectrum is to identify 
one peak as a Stokes process (RF creates a quasiparticle 
excitation) giving rise to the dissociation part of the RF 
spectrum and the other as an Anti-Stokes process (RF 
destroys a quasiparticle excitation). 

As the local imbalance is further increased beyond the 
superfluid to normal (SF-N) transition [22], see fig. [l]i, 
the majority spectrum no longer shows a local double- 
peak structure. This is consistent with theoretical work 
[23| [24] attributing the double peak structure in the nor- 
mal phase in previously reported RF spectra [25, 26 to 
the inhomogeneous density distribution. For increasing 
spin polarization the majority and minority pairing peaks 
lose spectral overlap. We interpret the missing overlap as 
indication that the minority atoms are no longer bound 
in pairs, each of them interacting with more than one 
majority atom, a situation we refer to as polaronic bind- 
ing. We have seen [22] that on the BEC side of the Fes- 
hbach resonance the overlap between minority and ma- 
jority spectra does not depend strongly on the presence 
of excess fermions. This is expected in a molecular pic- 
ture, where the pairing character is independent of the 
presence of excess fermions. At unitarity, within our ex- 
perimental resolution, the overlap starts to decrease at 
the superfluid-to-normal interface, see fig. This raises 
the question whether at unitarity a local probe of the mi- 
croscopic pairing, like RF spectroscopy, can distinguish 
whether the system is superfluid or not. 

Even when the spectral overlap decreases, there is still 
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FIG. 2: (color online) Spatially resolved RF spectra of an 
imbalanced Fermi gas at unitarity. a) The right half shows 
the majority spectra as a function of position in the trap 
expressed in terms of the majority Fermi radius Rj, the left 
half displays the minority spectra. The superfluid to normal 
transition region is marked by the gray vertical lines. The 
local polarization ai oc is given by the dashed red line. The 
error bars are the standard deviation of the mean value. The 
arrows indicate the position of the four spectra shown in fig. 
^ The image is a bilinear interpolation of 2500 data points, 
each plotted data point in the image is the average of three 
measured data points. The spatial resolution of the image is 
0.045-i? T 
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FIG. 3: (color online) Creation and spectroscopy of quasipar- 
ticles . a) Population imbalance thermally generates quasipar- 
ticles even at low temperatures comparable to A — ji^. /ij is 
the chemical potential of the majority component, b) The RF 
spectrum consists of a quasiparticle peak at negative frequen- 
cies and the pair dissociation spectrum at positive frequencies 
(dotted line). On resonance, the Hartree contribution U acts 
as an effective attraction and hence shifts the entire spectrum 
into the positive direction. 



equal response to the RF excitation in the high frequency 
tails, see fig. [TJ^ and fig. [T]i. These tails correspond to 
large momentum components in the interparticle wave 
function and hence address the short range physics. We 
expect this part of the spectrum to be insensitive to 
changes in the binding at large distances. 

The direct comparison between majority and minor- 
ity spectra clarifies our previous experimental results on 
minority RF excitation spectra in the |1) — |2) mixture 
[26] . in which we concluded that there is strong pairing 
in the normal phase. However, the observed spectral gap 
in the normal phase should not be interpreted as a signa- 
ture of pairing but rather as strong pairing correlations 
in the form of a polaron as suggested in [27j [28j [29] . 
The change in pairing correlations is indiscernible in the 
minority spectrum alone, but shows up in the spectral 
overlap with the majority spectrum. 

We now turn to a quantitative analysis of the spec- 
tral peaks in the superfluid phase for small density im- 
balance, and to the determination of the superfluid gap. 
Earlier work [15] |25] tried to determine the gap from the 
onset of the pair dissociation spectrum. However, the 
RF spectrum is not only sensitive to final state interac- 
tions, it is also shifted by Hartree energies [T2j [30] , as we 
show here. Furthermore, RF spectroscopy can excite all 
fermions, even deep in the Fermi sea [15] . Therefore, the 
onset of the pair dissociation spectrum occurs for atoms 
with momentum k = and, in the BCS limit depends 
quadratically on the gap parameter (co t h = The 
excitation gap can be directly observed if quasiparticles 
near the dispersion minimum are selectively excited, as 
in tunnelling experiments. 

Our solution is to study not the ground state of a 
superfluid, but excited states where quasiparticles are 
present. In a simple BCS description, quasiparticles are 
in pure momentum states, but increase the total energy 
of the system because their momentum state is no longer 



available to the other particles for pairing. Consequently, 
in an excitation spectrum, quasiparticles appear at neg- 
ative frequencies relative to the bare atomic transition 
frequency. The lowest energy quasiparticle appears at 
frequency —A, see fig. [3] 

Final state interactions and Hartree terms can also cre- 
ate line shifts, and two peaks are needed for analysis, the 
dissociation peak and the quasiparticle peak in our case. 
In essence, it is the separation between the peaks in spec- 
tra like fig. [TJd, which allows us to determine A. 

Thermal population of quasiparticles requires a tem- 
perature on the order of the excitation gap A. At uni- 
tarity, this temperature can be estimated to be 95% of 
the critical temperature, away from the low temperature 
limit addressed in this letter. Indeed, in samples of equal 
population of the spin states we were not able to spectro- 
scopically resolve any local double peak structures [22] . 
This problem can be overcome by introducing density 
imbalance between the constituents: The Fermi pressure 
(chemical potential /if) of majority atoms forces a fi- 
nite quasiparticle occupation into the superfluid region 
already at very low temperature. This allows us to se- 
lectively populate quasiparticles at the minimum of the 
dispersion curve see fig. [3^. 
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FIG. 4: (color online) a) Normalized peak positions of pair- 
ing peaks and quasiparticle peak at unitarity as a function 
of position in the trap. The SF-N boundary is marked by 
the dashed vertical lines. The arrow indicates the limit of 
low quasiparticle population relevant for b). Majority: blue 
open squares (pairing peak) and solid black circles (quasipar- 
ticle peak), Minority: solid red triangles, b) Pairing peak and 
quasiparticle peak positions as a function of the local interac- 
tion strength l//ci?a in the limit of small local imbalance (see 
arrow in a). Pairing peak: Solid circles, quasiparticles peaks: 
Open circles. 



In figure [4^ the position of the peaks of the majority 
and minority spectra are plotted normalized by e^t as 
a function of position in the trap for the unitary limit 
[22] . The peak positions are proportional to the local 
Fermi energy inside the superfluid region within our ex- 
perimental resolution. In the region of superfluidity with 
finite polarization the spectra show local double peaks. 
The position of the two peaks in the limit of small po- 
larization is depicted in fig. [4]o for various interaction 
strengths. 
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TABLE I: Superfluid gap A, Hartree term U and final state 
interaction E final in terms of the Fermi energy cf^ for various 
interaction strengths 1/kFd- 



1/kFd 


A 


U 


E final 


-0.25 


0.22 


-0.22 


0.22 





0.44 


-0.43 


0.16 


0.38 


0.7 


-0.59 


0.14 


0.68 


0.99 


-0.87 


0.12 



It was unexpected that the quasiparticles appear at 
positive frequencies (relative to the atomic transition fre- 
quency) . This is caused by the presence of Hartree terms 
[12], resulting in an overall shift of the systems energy 
and the RF spectrum [22]. In the weakly interacting 
limits, the Hartree term reduces to a simple mean field 
shift. In the strongly interacting regime one has to resort 
to QMC calculations [30j [3TJ [32] for a numerical value of 
U. 

In a mean-field description of the balanced superfluid 
starting from the BCS-Leggett ansatz for the BEC-BCS 
crossover [18 taking into account the Hartree term J7, the 
dispersion relation of the quasiparticles can be expressed 

as E k = ^A 2 + (e fc + /7-/i) 2 [T2] where e k = ^ is 
free particle kinetic energy and \i is the chemical poten- 
tial. This mean-field formalism gives the analytic expres- 
sion for the two peak positions. A quasiparticle at the 
minimum of the dispersion curve will respond at an RF 
offset of uj rf = -E kmin - fi + e kmin = -A - U, and 
the maximum of the pair dissociation spectrum occurs 

at huu max = | (yV 2 + ^A 2 -//) - U ~ ±<Jth - U, 
where // = \i — U and uo t h is the dissociation threshold 
(which is at momentum k = 0). 

We determined the superfluid gap A and the Hartree 
energy U from the peak positions in the limit of small 
density imbalance (cri oc ~ 0.03). At unitarity with the 
chemical potential \i = 0.42 e^, confirmed in previ- 
ous experiments, see [33] and references therein, we ob- 
tained A = 0.44 cf^ and U = — 0. 43 epj, in excellent 
agreement with the predicted values A t — 0.45 and 
Ut = — 0.43ei?t from QMC calculations [34]. Our de- 
termined values for A and U values suggest the mini- 
mum of the quasiparticle dispersion curve to occur at 
k m i n ~ 0.9kp- Table [I] shows the gap and Hartree energy 
for various interaction strengths. Away from unitarity 
we relied on QMC calculations for the chemical potential 

At ESI- 

For an accurate quantitative comparison [22 final state 
interactions, also listed in table [T| had to be taken into 
account. The effect of final state interactions is an overall 
mean field shift of Ef ina i = 47T ^ a n. This shift affects 
both the quasiparticle peak and the pairing peak equally. 



In conclusion, we have performed spatially resolved RF 
spectroscopy of majority and minority components of a 
trapped imbalanced Fermi gas in the strongly interact- 
ing regime with small final state effects. In crossing the 
superfluid to normal boundary we observed a gradual 
crossover in the pairing mechanism by comparing major- 
ity and minority spectra. The majority spectrum shows 
a local double peak spectrum in the polarized superfluid 
region which allowed us to determine the superfluid gap 
A and the Hartree terms U. The spectra in the normal 
phase are consistent with a polaron picture. 
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ing discussions and A. Sommer for critical reading of the 
manuscript. This work was supported by the NSF and 
ONR, through a MURI program, and under ARO Award 
W911NF-07-1-0493 with funds from the DARPA OLE 
program. 



Auxiliary Material: Determination of the Superfluid 
Gap in Atomic Fermi Gases by Quasiparticle 
Spectroscopy 

Determination of the superfluid boundary 

It has been shown previously [8 that at unitarity the 
difference column density profiles serve as an indicator 
for the SF-N boundary. The discontinuity of the minority 
density results in a pronounced "cusp" in the difference 
profile n^(r) — n^(r), see fig. [5] In the main body of the 
paper the phase boundary has been determined by this 
peak position. Hence, "normal" refers to spatial regions 
beyond the peak, "superfluid" refers to spatial regions 
inside. 

The color coding in the graphs in fig. [5]shows where the 
spectral overlap (definition see below) between the ma- 
jority and minority pairing peaks is lost. On resonance 
(B = 690 G) this position shows excellent agreement with 
the position of the cusp in the difference density profile 
and can therefore serve as an alternative indicator for 
the SF-N boundary. This coincidence breaks down away 
from resonance: On the BCS side (B = 710 G) the spec- 
tra are less "robust" against polarization and spectral 
overlap is lost before the column density difference shows 
a peak. The reverse situation occurs on the BEC side of 
the resonance (B = 671 G). Note that on the BEC side 
the minority cloud does not extend much further than 
the peak position in the column density difference. 

Quantification of spectral overlap 

In order to calculate the spectral overlap of a spec- 
trum, the quasiparticle peak was fitted by a gaussian and 
subtracted from the majority spectrum. The overlap is 
then defined as one minus the difference of the integrated 
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FIG. 5: (Color online) The difference column density profile 7ij(r) — nj,(r) (radially averaged), a) BCS side (B = 710C7), b) 
Unitary limit (B = 690C7), c) BEC side (B = 671C7). Each profile shown is the average of 10 individual profiles. Blue marks 
the region of complete spectral overlap between majority and minority components, red marks the region where there is no 
complete spectral overlap. The black arrows at the bottom indicate the radial size of the minority component. 
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FIG. 6: (color online) Overlap of majority and minority pair- 
ing peak as a function of position in the trap for various in- 
teraction strengths. A power law was fitted to the curves as 
a guide to the eye. Unitary limit (black circles): l//ci?t a = 0> 
phase boundary at r c /R^ — 0.46; BEC side (red triangles): 
1/kpia = 0.39(1), r c /i? T ~ 0.45; and BCS side (blue dia- 
monds): l/kp^a = -0.25(1), r c /i? T ~ 0.35 



spectra normalized by the sum of the integrated spectra, 
see fig. [6j As mentioned in the main text and above, on 
the BEC side of the Feshbach resonance almost complete 
spectral overlap can be observed into the normal region. 
On the BCS side the reverse situation occurs. 

Minority peak shift for high imbalance 

Fig. 4 in the main body of the text shows the peak 
positions normalized by the local majority Fermi energy 
e^t- Fig. [7] shows the bare peak positions as observed in 
the experiment. 

One unexpected finding in fig. 4a in the main body of 
the text is the sudden increase of the minority peak po- 




FIG. 7: (color online) a) Local majority Fermi energy in kHZ, 
the error bars are the standard deviation of the mean value, 
b) Peak positions of majority and minority in kHz. Major- 
ity (Blue): Pairing peak (higher frequencies, only discernible 
in the SF region) and quasiparticle peak (lower frequencies). 
Minority (red): The peak can be traced well into the normal 
region. 



sition for high imbalances. This behavior can be traced 
back to the data in fig. [7| The minority peak position 
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(red) shows a change in curvature as the SF-N bound- 
ary is crossed, in contrast to the majority Fermi energy. 
Therefore, the ratio of minority peak and majority Fermi 
energy shows a sudden increase towards the edges of the 
cloud. The value of the peak position of uo po i c± 0.9 ep^ 
(where final state interactions of Efi na \ ~ 0.05 ep"\ have 
been taken into account) is higher than the theoretically 
calculated value of ~ 0.6 ep^ 0. 



Quasiparticles in an equal mixture 

In a previous publication [16] low temperature RF 
spectra did not show any signatures of quasiparticles. We 
have attempted to create thermally generated quasiparti- 
cles in an equal density mixture for higher temperatures 
T/Tp > 0.20. The experimental results in fig. [8] show a 
decrease of the gap parameter, but no local double peak 
structures could be resolved for any temperature. 

We can estimate the temperature required to populate 
quasiparticles, assuming that the temperature has to be 
on the order of the gap and that A(T) is given by the BCS 



T r ? ]. At unitarity, A /T c 



relation = 1.74y 1 T 
3. One would expect to populate quasiparticles only very 
close to the transition temperature T/T c ~ 0.95, when 
the gap is only one third of its low-temperature value. 



The theoretical dissociation spectrum including the 
Hartree term U 



Starting from a BCS-Leggett mean-field wavefunction 
and applying Fermi's Golden Rule, the RF spectrum can 
be described by (neglecting the Hartree term): 



Resolution / Experimental broadening 

For comparison with the theoretical spectrum, we 
model the RF pulse of T = 200 fis length as a square 
pulse, in the frequency domain, resulting in a FWHM of 
the RF spectral power of Av = 2 • ^ 1.39 f ~ 4.4 kHz. 
The theoretical spectrum consists of two parts: 1) A dis- 
sociation peak including the Hartree energy, described by 
eqn. [I] with A and U as given in table I in the main text. 
2) A quasiparticle peak modeled as a narrow (FWHM 
= 1 kHz) Lorentzian with a peak height adjusted so that 
it resembles our data. This spectrum was convolved with 

the Fourier transform of a square pulse f(uo) oc 



m 2 



Fig. [9] shows that the theoretical spectrum reproduces 
our data quite well. The deviation in [9}d) might be at- 
tributed to additional broadening mechanisms like finite 
quasiparticle lifetime, finite temperature and atomic dif- 
fusion during the duration of the RF pulse. The convolu- 
tion causes a small shift of 0.05 ep-\ in the spectral peak 
position due to the asymmetry of the theoretical spec- 
trum and has been accounted for in the determination of 
A and U. 

The value given above for the experimental resolution 
is confirmed by looking at the blurring of the sharp onset 
of pair dissociation. Equation [l] predicts that the thresh- 
old and peak position in the strongly interacting regime 
differ by less than 10% of the Fermi energy. Adjusting 
the experimental resolution to ~ 4 kHz accounts for the 
experimentally observed difference of threshold and peak 
position of about 0.3 epp 



r (u oc — wi + — + — 1 

uo z V uo uo 



where uj t h = \/ A 2 + [i 2 — fi is the dissociation thresh- 
old and \i is the chemical potential. 

The corresponding quasiparticle dispersion relation is: 



(2) 



E k = Va 2 + (e, - n)* 

e k being the free particle dispersion e k = 

Hartree terms U modify the quasiparticle excitation 
spectrum [151150]: 



E k = y^A 2 + (e k - (/i - U)) 2 
resulting in an RF spectrum of 



LO u V CO 



2n> 



(3) 



(4) 



where a/ = uo + J7, uu' th = Uth + U and \jl = /i — U. 
This demonstrates that the spectrum retains its func- 
tional form but the entire spectrum is shifted by U. 
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FIG. 8: (Color online) Local RF spectra of an equal spin mixture for various normalized local temperatures T/Tpf. a) 
T/Tp| ~ 0.20, b) T/Tp| ~ 0.22, c) T/T FT ~ 0.34, d) T/T FT ~ 0.55. No local double peak spectrum can be resolved in the RF 
spectrum. For comparison, the double peak spectrum of an imbalanced mixture with T/Tp^ — 0.06 is added with a red dashed 
line. 
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FIG. 9: (color online) Comparison of experimental (red) and 
theoretical (black) line shapes for spectra a) and b) of figure 
1 of the main body of the paper. The theoretical curve is 
obtained from a BCS-Leggett mean field description including 
the Hartree term and a convolution with the experimental 
resolution of ~ 4.4 kHz. The values for A and U as calculated 
from the peak positions lead to a reasonable agreement with 
the data. 
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